ON THE SOLUTION OF LINEAR DIFFERENTIAL EQUA¬ 
TIONS OF SUCCESSIVE APPROXIMATIONS. 


By PRESTON A. LAMBERT. 

(Read April 20 , 1911 .) 

The object of this paper is to apply to the solution of linear differ¬ 
ential equations, both ordinary and partial, the method of expansion 
into series used in the solution of algebraic equations in the papers 
read bv the author before the Philosophical Society in April, 1903, 
and in April, 1908. 


Let the given differential equation be 


(1) 


/ 


( . , d y_ dl y a "y \ _ 

\ ’ dx ’ dx"-' ' dx" ) 


o. 


The method of solution consists of the following steps: 

(a) Break up the left-hand member of the differential equation 
into two parts. 


and 


/, 


A 


dy 

d-y 

d"y 

dx ’ 

' dx-' ' 

'' ’ dx" 

dy 

dy 

d"y 

dx 1 

' dx 1 ' ’ 

' ‘ ’ dx" 


such that the first part equated to zero can be integrated by some 
known method, and multiply the second part by a parameter S, inde¬ 
pendent of .1* and y. Replace the given equation by 


, . f ( dy dy d“y\ ( dy dy d“y\_ 

( 2 ) /,(*,/> ({x , dx i, • • •. dx 2 ) + S ^\ X ' y ' dx' dx 1 ' " dx") 

( b) Assume that 

(3) y 0 + \\S + y. 2 S 2 + yS A -)- y 4 N 4 + • • • 
makes equation (2) an identity. 


O. 
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(c) In this identity arranged according to the ascending powers 
of S equate to zero the coefficients of the different powers of S. 

( d ) Solve the differential equations thus obtained in regular order 
for y 0 , 3T3T JT 3T 

(c) Substitute these values in (3 ) and make S unity. The result¬ 
ing value of v, if it contains a finite number of terms or if it is a 
uniformly convergent infinite series, is a solution of the given dif¬ 
ferential equation. 1 

The method of solution of linear differential equations as here 
outlined does not seem to occur in mathematical literature except as 
developed by the author. 

The method will be exemplified by applying it to two differential 
equations, important in mathematical physics—Bessel’s equation, a 
second order ordinary differential equation, and Fourier's equation 
for the flow of heat, a second order partial differential equation. 

Bessel’s equation is 



Replace Bessel's equation by 



and assume that 


y — Vo + y\S + v 2 ^ _ + yzS 3 + + • • • 


makes the latter equation an identity. 

When arranged in ascending powers of 3 T this identity is 



d~j '2 S' T • • • = o. 


1 This method gives a formal solution of non-linear differential equations, 
but up to the present time the author has been unable to test the resulting 
series for convergency. 
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Equating to zero the coefficients of the powers of 5 * in this identity, 
there result the following differential equations for the determina¬ 
tion of y 09 y 19 y 2 ,y i9 


„ 2 d\y n dr n 

' l dx 2 +X dx~ n ^ = °’ 

” 2 - /j 1 + t - >?y x + *y 0 = °- 


dx l ^ dx 
2 , d > 


.,1 “ S2 


+ x-~— try 2 + xy x = o. 


dx 2 ' dx 


The equation in y 0 is a homogeneous linear differential equation 
and its solution is 


y 0 = Ax n + Bx~ n . 


Substituting this value of y 0 the equation for determining y x 
becomes 


d\)\ 

dx 2 


+ -t' 


dx 


n 2 j\ = — /Er' l+2 — Bx 71+2 . 


This equation becomes exact when multiplied by x~ n ~ l . The 
resulting equation integrated gives a linear equation of the first order, 
the solution of which is 

— Ax n+2 Bx ~ n+2 

^ l ~ 2 2 (;/+ i) 2 2 (/* — l) 

Substituting this value of \\ in the equation for determining y 2 
and proceeding in the same manner 


y 0 = 


2 4 


Ax ”+ 4 

2 ! (tl + I ){n + 2) 2 4 • 2 ! 


B, r~ n+4 

(n _ !)(;/ - 2) 


In like manner 


- Ax n+6 


y$ 


2 6 *3 !(« ’+ 1 )(« + 2 )(« + 3 ) + 2 6 -3 !(« - 0 (« — 2 )(« - 3 )’ 


B.r ~ n+6 


and so on. 
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Substituting these values of 3V 3V V 2 , y 3 , •••in 

y = Vo ~h y\S “I - 3 ^ + v 3 ^ 3 ~h V-^ 4 + 

and making 6* unity, 

r i x 2 i _ & 

* X [_* n -f i 2 2 (;/ + i)(y/ + 2) 2 4 • 2! 


(;/ + i)(// + 2)(;/ + 3) 2 6 .3 ! 

T i .r 2 i .r 4 

“h Bx M ] I T* o “1" / w r —4 . 

|_ // — I 2 2 (;i — 1)(// — 2) 2 4 • 2! 




+ ■ 


x° 


+ •• 


] 


(« - 0(« - 2 )0 — 3) 2 6 ' 3 ! 

When n is not an integer the terms of both series in this value 
of 3' continue indefinitely according to the law of formation which 
inspection makes evident, both series are uniformly convergent ex¬ 
cept when .r = o, and both series are solutions of the given differ¬ 
ential equation. 

When 11 is a negative integer the law of formation of the terms 
of the first series changes after the (;/)th term and when n is a 
positive integer the law of formation of the terms of the second 
series changes after the (;j)th term. The second case will be con¬ 
sidered. 

When n is a positive integer the (;/)th term of the second series is 

Bx n ~ 2 

fn-i ~ 2 a»— I)(„ _ 1)! • 

Substituting this value of y„_i in the differential equation for 
determining y n , 

. 7 d 2 y dy „ 

and solving for y„ by the method used in solving for 3V3V3V “*> 

a= 2 -. ; 7rf« - t): L' 1 '" 

I11 determining y,i + i. v„_ 2 , y„ i3 , •••, the second term in the bracket 
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gives the terms of the first series in the value of y multiplied by a 
constant. This new series is combined with the first series in the 
value of y. 

The first term in the bracket gives 

- B I" x"+ 2 log x x n+2 / i \ "I 

— 2 2n ~ l n ! (;/ — i) ! [_ 2 \)i + i) 2 2 {n + i) \ 1 + n + i / J ’ 

- B r x"+ 4 log x 
Jn+2 — 2 u ~ l n ! (;/ — i)! [ 2 ! 2 \?i + i)(;/ + 2) 

x’ t+1 /11 1 \1 

~ 2 ! 2 \)t + 1 )(n + 2) \ 1 + 2 n + i + n +2) y 


The solution of Bessel's differential equation when n is a positive 
integer is therefore 


'= Ax n r 1 — 1 

L n + 1 2“ ~ 


I X 

(;/ 4- i)(« 4- 2) 2 4 -2 ! 

1 


0 + i)(a + 2)(« + 3) 2 6 *3 ^ 


+ 




+ Bx~" I 1 + 1 . 1 2 + 


// — I 2 2 (U — l)(^ — 2) 2 4 - 2 ! 


_ x°~ n ~~ I 

+ + 2 2(n ~ l \u — i)!(//— i)! J 


/)Viogj r 1 x 2 
~~ 2 2 " -1 ;/1 (n — 1)! L* n + 1 32 


1 x 4 

(;/ 4- 1)(// -f 2) 2 4 • 2 ! 


1 


+ 


(u + i)(;/ + 2)(« + 3) 2 6 - 3 ! 

/xi' n r 1 ( 1 \ ^ 

2 2 ,i_ 1 7 / ! (n — 0 1 L ;/ + 1 \ ^ + 1/ 2 2 

I / I I I \ X 4 “I 

“ (n + 1)(// + 2) \ 1 2 + « + 1 + ;/ 4- 2 / 2 4 * 2 ! + ‘ ‘ ‘ J * 

This is also the solution of the differential equation when n is a 
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Fourier's partial differential equation for the linear flow of 
heat is 

dV d-V 

W~ A eP' 

Replace Fourier's equation by 

dV Q d 2 V 

dt ~ Sk dx 2 

and assume that 

v=v 0 + r 1 s+v 2 s'- + v 3 s° + --- 

makes the latter equation an identity. 

When arranged in ascending powers of S this identity is 

! S 3 + ... = o. 


d K\ 

dt 

dv t 
+ ~dT 

s + W 

S ,, dV -i 

5 + ~dF 


7 dx- 

dx- 

A dx 2 


Equating to zero the coefficient of the powers of S in this identity, 
there result the following partial differential equations for the deter¬ 
mination of V 0 , J\, Vo, V z , • * 


S K 


-o. 


dV >_ K dl v x 


dV *_K d V^-o 

dt A dx 2 ~ ’ 


dt 
J-V 

dt ~ “ dx 1 

These partial differential equations solved in regular order give 

k - K=r(-m, k~ f. 

Substituting these values of V 0 , V lf V 2 , ... in the assumed 

value of V and finally making S' unity, there results 

(A) V- <f>(. r) + 4 >"(x)(Kt) + 4 > lv (.v) ® + </>'» 
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which is a solution of Fourier’s equation for all values of </>(.r) for 
which V either contains a finite number of terms or is an infinite 
series uniformly convergent both in .r and in t . 

The following table shows several values of <f>(x) and the cor¬ 
responding solutions of Fourier’s equation. 


I (i) <p(x) = A, 

(2) 4 >(x) = Ax, 


V=A, 
V= Ax, 


(3) = Ax\ 

(4) <f>(x) = A sin (fix), 

(5) </>(F) = A cos (nx), 

(6) 4 >(x) = Ae"*, 

(7) <X- r ) = Ae ~ n \ 


V= A(x 2 + 2 A?), 

V = Ae~ n ~ Kt sin {nx) y 

V = Ae~ n * Kt cos 
V= Ae" x +"* Kt , 

V = Ae- nx+n2Kt , 


(8) 4>(x) = Ae nx sin (nx), V = Ae nx sin (iix -fi 2 /rKt). 

It will be noticed that in these solutions <j>(. r) is the value of V 
when t = o, that is V = <f>(x) is the initial heat distribution. 

It will also be noticed that in all these results x may be replaced 
by xa. This statement is true of the results in the several fol¬ 
lowing tables. 


If Fourier’s differential equation is replaced by 

BV_ d*V 

S dt - K dx- 

and the assumption made that 


v=v 0 + r 1 s’ + r 2 s ’ 3 + r 3 s*+..- 


makes this equation an identity, this identity arranged in ascending 
powers of .S' is 

K dx 2 


cPV, 

d 2 K 

d 2 K 1 

+ K d J 

5 + A a? ? 



dV, 

dK 

dt 

dt 

dt 


Equating to zero the coefficients of the powers of .S' in this 
identity, 
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d 2 F. d- V. I dV. 8 2 V„ 1 8 V. 

8 hF 8 x 2 ~ K ~ 8 t ~ °’ 8 x 2 ~K Tt = °’ " 

Solving these partial differential equations in regular order for 
V o, V 2 , V z , •••, substituting these values in the assumed expres¬ 

sion for V, and finally making S' unity, the result 

r= 4>(t)x + g <f>'(t) -j + ^- 2 4>"{t) -.-j + • • • 

(R) ’ ’ 4 

+ 6(0 + y e '(0 2 , + b 6 "^. + ’ ’ ‘ 

is a solution of Fourier's differential equation for all values of 
and 6 {t) for which V either contains a finite number of terms 
or is an infinite series uniformly convergent both for x and for t. 

Solutions of the differential equation when </>(t) =0 correspond¬ 
ing to several values of 0(t) are as follows— 


II 4>(t) = o, 

(o m = 

( 2 ) e{t) = At, 

(3) 0(t) = At 2 , 

(4) 8{t) = AR, 

( 5 ) 0(t) = Art, 

(6) 6(t) = Afi, 


( 7 ) 6(f) = Art, 


V=A, 


V=A { ( + Jjc)> 

r-A^ + F + Ttf*)- 

v= ar[ 


x- x 

2 !(2A 7 ) ~ 4 \( 2 Ktf 

t -6 

+ 


_ 3- r 
6 ! (2 A?) 3 


V= 


A e ia 

KR ’ 



V=At 


t 


3 - 1 ' 2 3 x> 

2 ! 2A? + 4I (2A?) 2 



3 

2 ! 


3 

6 ! (2 A?) 3 



x 2 3*5 x 4 

zKt 4 ! (2AV) 2 


3*5*7 

— (W + 



> 
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I U X 

(8) 9 (t) sin (///), V — A\ sin (;//) -f ^ cos nt -y 

n 2 -r 4 *1 

- K , sm (nt) 4 ,-J, 

( 9 ) = Alogt, V = A ^log t + ~ A- 

2 X 6 2 • 3 X s 1 

+ K 3 t 3 61 ~ A" 4 / 4 8! + 'J’ 

CO) m-Ar\ V-AS-[ , + * ' + £ + ...]. 

It will be noticed that in these solutions V = 6 (t) is the heat dis¬ 

tribution when x = o. 

Solutions of the differential equation when 0 (t) — o correspond¬ 
ing to several values of </>(/) are as follows: 


III 0(/) = O, 

(1) 4>(f) = A, V= Ax, 

(2) = At, J = A j^.t/ + j j , 

T „ 2.V 3 / 2.1' 5 n 

( 3 ) = At~ l = A j^.r/ -f ^ > 

(4) ^ |^' r + ^ | — 2 2 A' 2 / 2 "JT 

3 - r? 

+ 2 3 A ^ 7 , ~ • 


(5) </>« = At~K 


(6) <£(/) = Ac" 1 , 

(7) ^>(/)=^sin(«/), 


F= 

At 

-*h 

1 

2 Kt 

-t' 3 3 

3 ! + ~?TOP " 





3 * 3 
2 3 A"¥ 

X 7 

7! + 

F= 

Ac 1 

r 

n w 3 

ir x rj 


■' 1 -V + 

K 3! 

+ A- S ! 

+ • • • 


r 


n 

cos (;//) 

.V 3 

V= 


sin (;//) 

V+ K 

3! 



1911-1 


LINEAR DIFFERENTIAL EQUATIONS. 


283 


(8) </>(/) = A log t, V= A |x log / + - j , 

I X 5 2 X 7 1 

~lCA "5] + 7i ~ ' ’ J 

It will be noticed that in this set of solutions T = o when x = o. 


Let u 1 = f 1 (.r, t), u 2 = f 2 (y,t), u s = f s (s,t) represent solutions 
of the three one-dimensional Fourier’s equations, 

dv_ , d 2 V dV_ a 2 _F dV d-v 

~dt = A d 7 ~’ ~dt K ~df ’ ~dt ~ ' dz 2 
respectively. It is readily proved that 

J r —u l u 2 and V = u^uit 3 


are solutions respectively of the two-dimensional Fourier’s equation 


dV T Jd 2 V 
dt ~ \ dx 2 



) 


and the three-dimensional Fourier’s equation 


dV _ fd 2 V d 2 V d 2 V\ 
dt — dx 2 dy 2 dz 2 J 

This shows how solutions of the two- and three-dimensional 
Fourier’s equations can be obtained from the solutions of the one¬ 
dimensional equation. 

For example, from the one-dimensional solutions 

X- 

zJ p 4Kt 

V = — yx ~\ —■ and V — Ae~ n ~ Kt sin (nx) 


the three-dimensional solutions 


IV (1) 


v = 


Ac 


K"'A 


(2) V = /!c- (a ' + P"~ l " r ' jK ' sin (a.r) sin (J 3 y) sin (7-), 
respectively, are obtained. 
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If the solution of the three-dimensional Fourier's equation 

dV / d 2 V FT d 2 F\ 

~W = K \~d7 2 + ~df + ~d?~ ) 

is a function of r and t only, so that 

V = f(r,t ), where r = (a* 2 -|- v 2 + ^ 2 )-, 

the transformation of the given equation from rectangular to polar 
coordinates shows that the solution is 



where u is a solution of the Fourier's equation 

du ^ a 2 ;/ 

5 t ~ A dr ' 

It follows that solutions of the three-dimensional equation of the 
form V = f(r , t ) are obtained by replacing .r by r in any solution 
of the one-dimensional equation 

dV _ K dht 
W “ A dx 1 ’ 

and dividing the result by r. 

In this manner are obtained the solutions 



Ae 4Kt 
Tl<7* ’ 

A 

(3) I w = — e nr sin {nr -f- 2 irKt). 

It is interesting to compare the solutions of Fourier's partial dif¬ 
ferential equation obtained in this paper with the solutions tabulated 
by Sir William Thomson in the mathematcal appendix of the article 
on “ Heat " in the Encyclopaedia Britannica,” ninth edition. 

Sir William Thomson obtains his results by summation, that is 


(2) v= 
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by integration, from the solution IV (i) above. All his results 
occur directly in the above tables or are combinations of two of these 
solutions. It is evident that there are several misprints in the results 
as printed in the “ Britannica.” 

Of course there are many solutions of Fourier’s equation which 
must be built up from elementary solutions, however found, by 
means of Fourier series, or which must be obtained by the methods 
of harmonic analysis. 

The solution III (5) above is the series used by Sir William 
Thomson in his solution of the problem of the secular cooling of the 
earth. 2 

An interesting result in pure mathematics is obtained as follows: 
Sir William Thomson shows that for a continued point source of 
heat, if the rate is an arbitrary function of the time, /(f), the solu¬ 
tion of Fourier’s equation when K = 1 is given by the definite 
integral 



The second part of the general solution (B) above shows that 

is also the solution of Fourier's equation for the same conditions. 

It follows that 


r™ - ■'"> - 4b[7 m TT + + ] 


is a general formula for computing the definite integral. 

Lehigh University, 

Bethlehem, Pa. 

2 “ Mathematical and Physical Papers/' Vol. III. 


